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Abstract 

We propose a simple scaling theory describing the variation of the mean first passage time (MFPT) t(N,M) of a 
regular block copolymer of chain length N and block size M which is dragged through a selective liquid-liquid interface 
by an external field B. The theory predicts a non-Arrhenian r vs. B relationship which depends strongly on the size 
of the blocks, M, and rather weakly on the total polymer length, TV. The overall behavior is strongly influenced by the 
degree of selectivity between the two solvents \. 

The variation of t(N, M) with TV and M in the regimes of weak and strong selectivity of the interface is also studied 
by means of computer simulations using a dynamic Monte Carlo coarse-grained model. Good qualitative agreement with 
theoretical predictions is found. The MFPT distribution is found to be well described by a P - distribution. Transition 
dynamics of ring- and telechelic polymers is also examined and compared to that of the linear chains. 

The strong sensitivity of the "capture" time t(N, M) with respect to block length M suggests a possible application 
as a new type of chromatography designed to separate and purify complex mixtures with different block sizes of the 
individual macromolecules. 

1 Introduction 

In a recent series of studies [1-4], we have reported the results of comprehensive computer simulations which model 
the behavior of hydrophobic - polar (HP) multiblock copolymers at a selective penetrable interface which divides two 
immiscible liquids, like water and oil, each of them being favored by one of the two types of monomers. As is well known 
from experiment [5-8] , in the presence of selective interfaces for which the energy gain for a monomer in the favored solvent 
is large, the hydrophobic and polar blocks of a copolymer chain try to stay on different sides of the boundary between the 
two solvents, leading thus to a major reduction of the interfacial tension which has important technological applications, 
i.e. for compatilizers, thickeners or emulsifiers. Not surprisingly, during the last two decades the problem has gained a lot 
of attention also from theory [9-13] as well as from computer experiment [14-18]. While in earlier studies attention has 
been mostly focused on diblock copolymers [6,8] due to their relatively simple structure, the scientific interest shifted later 



to random HP-copolymers at penetrable interfaces [10-13,18]. In contrast, our investigations have focused mainly on the 
impact of block size M on the static properties and localization kinetics of regular multiblock copolymers at the phase 
boundary between the two immiscible solvents. We showed that these are well described by a simple scaling theory [1,3] 
in terms of the total copolymer length N (the number of repeating units in the chain) , the block size M (the number of 
consecutive monomers of the same kind), as well as the selectivity parameter that is, the energy gained by a monomer 
when moving into the more favorable solvent. 

In the present work we extend our investigations to field-driven multiblock copolymers which cross the liquid-liquid 
interface and are temporarily trapped for a typical time t(N,M). Our simulation results reveal a rich behavior of the 
drifting chains with qualitative differences in the two regimes of weak and strong selectivity x- We demonstrate that a 
scaling theory may be designed which captures the main qualitative features of the Mean First Passage Time t(N, M) in 
terms of chain length N and block size M and provides insight into the complex behavior of the copolymer being dragged 
through the interface. We will show that the time t(N, M) is very sensitive towards the block length M. This results can 
be viewed as a potential ground for the possible development of block - sensitive chramotography designed to separate 
complex copolymer mixtures. It is worthy of note that the dynamic Monte Carlo simulation, used in our study, does not 
take into account the hydrodynamic interactions. The role of hydrodynamic effects might be very important and deserves 
future investigation. 

In SectionElwe introduce our model and derive the main analytical results which describe the variation of translocation 
time with block size M and force strength B within the framework of Kramers approximation for the first passage time. 
The main features of the Monte Carlo model which is used in our simulations are briefly sketched in Section 13 . II while in 
Section [3.21 we discuss the most salient results of our computer experiments. We end the paper with the main conclusions 
which the present investigation suggests. 

2 Theory of copolymer detachment from the interface 

We now consider the model of translocation dynamics where a growing fraction of the chain is detached from the interface 
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Figure 1: Schematic representation of the string of blobs subject to an external force B. In this scheme, the leading part 
of the chain (dangling end) is growing at the expense of the string of blobs still captured at the interface. P and H denote 
the polar and hydrophobic scmispaces respectively. 
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by the external field B and dragged into the underlying solvent as shown schematically in Figure ^ The physical picture 
which provides the basis for our detachment model will be justified in Sec. 13. 21 by the direct visualization within the Monte 
- Carlo simulation. This model will be used here to describe the characteristic time r it takes to tear a string of blobs off 
the interface. To this end we calculate the free energy function. 

2.1 Free energy of the detachment model 

We consider the free energy of the model, depicted in Figure^ as a function of the instant length of the chain, m, which is 
in the lower half (H) of the container (the dangling end length). The total free energy then contains the following terms: 

F{m)= F scl (m) + F cl (m) +F str i„ g (m) . (1) 
v * ' 

dangling end free energy 

In this equation F se i(m) and F e \(m) are the selectivity interaction and the elastic deformation contributions of the leading 
end, whereas F s t r i ng (m) is associated with the energy of the remaining string of blobs. 

The selectivity energy of the dangling end is related to the interface potential and the HP - copolymer binary sequence 
<t(s) (c(s) assumes the values +1 and —1 depending on the nature of the monomer, see ref. [1]. In the hydrophobic 
semispace H (z < 0) this free energy reads 

N 

F se i(m) = ^ J dsa(s) (2) 

N—m 

A simple binary sequence for a multiblock copolymer with a block length of M can be constructed as <j(s) = sin (-^ s). 
Substituting this in eq|21yields: 

TP I \ * M ■ 2 ( 71171 \ fo\ 

Fsel{m ) = -— sm [—) (3) 
where we have also taken into account that N and M have been chosen as multiples of 2 (as is usual in MC - simulations). 

2.1.1 Free energy of elastic deformation 

Here we consider the deformation energy of the different parts of the chain based on scaling arguments. We assume the 
interface position to be at z — 0. The external field acts uniformly on each polymeric segment with a force B which points 
in the negative z direction. 

Each instantaneous configuration can be thought of as an array of loops which spread on both sides of the interface. 
We denote them as positive (z > 0) and negative (z < 0) loops respectively. The two types of loops respond differently 
to the external field B. The positive loops are compressed whereas the negative loops are stretched with respect to their 
equilibrium shape for B = at the interface. Consider first the elastic free energy due to compression of a positive loop. 
Suppose a loop of a length n runs in the positive direction up to a height h. Then the elastic free energy can be written as 

h 
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The first term in the r.h.s. of eq^Jis due to confinement entropy whereas the second one is related to the external field 
energy, T stands for the temperature and a denotes the Kuhn segments length. Optimizing eq0]with respect to h, one 
gets ^ = ^J!^Y which substituted back into eq^Jgives 

F«(n)~T^J n 2,/(i +2 ,) (5) 



Let us now look at the stretching free energy associated with the deformation of the a negative loops. The elongation 
of a chain in a good solvent due to forces / applied at both ends [19], for small extensions, an" < h <C an, reads 

a 2 n 2 »f/T , f^T/an" 

an{af /T) {1/v) - 1 , / > T/an v 

Following F. Brochard - Wyart et al. [20-22] (see also the Monte - Carlo simulation in the wake of these results [23]), one 

dh(s) ( f \ 

may describe the non-uniform stretching of tethered chains in the differential form, ~ a I 4? ) , where s is 

the running curvilinear coordinate of a polymeric segment and the field B acts on each monomer. The total force is then 
/ = sB. The elongation of a string with n monomers is then obtained from dh ^ — a (-^p 2 -)^"' 1 by integration of s, 
< s < n, and yields 

h * an* ( "4) ■ (7) 



, T J 

This result is valid for fairly large forces nB > T/an u , or B > T/n v+l , which exceed the thermal agitation. On the other 
hand, the Pincus' blob size can not be smaller than the Kuhn length, i.e. nB < T/a, or B < T /an. As a result, the 
deformation law JJJ is valid for forces in the interval T/an u+1 < B < T/an. The corresponding free energy of stretching 
takes the form 

i^")(n)~ J/f ds-a^f'^ 1 b] s^Vals-T^) 1 '"^ 1 '^ 1 . (8) 
o o 

In the weak field regime, B < T /an u+1 the elongation law follows from the first line of eqEl i.e. dh j^ ~ a (^jr) s 2v 
and the resulting stretching free energy becomes 

n n 

F(")(n)~ f ff s ds~a(^)Bjs' 2v + l ds~T{^) 2 rf v + 2 . (9) 



Thus, the the elastic free energy in both weak, and relatively strong stretching regimes (see eqs [5] and [S] ) , is given by 
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2.1.2 Total free energy 



>From the previous section we have for the elastic free energy of the leading part which includes m monomers of the 
dragged chain, according to eaUUI 

FexM^-T^y^m^Hi . (11) 

The negative sign here implies that the chain moves downhill along the external field. 

Consider now the free energy of the string of blobs, that is, of the rest of the copolymer which is still trapped at the 
interface. The length of such a blob at the interface is [1] 

g ^ x -V(l-u) M -{l+u)/(l-u) ^ (12) 



whereby the blobs are placed alternatively in the P - or H - sides of the interface. Each blob has an energy of the order of 
T. The number of blobs is (N — m)/g with g - the number of (both H and P) monomers in a blob. Taking into account 



the expressions for the elastic energy (JjjJ and i|l(J|) , the free energy of the string of blobs is proportional to 
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where ai and a.2 are some numeric factors and 7 = 1/(1 + 2z/). The first term in ea 1131 is associated with the energy of 
localization (therefore it is negative) whereas the two other terms are due to the elastic deformation of blobs. 

In the strong localization limit the positive and negative loops which constitute the string of blobs are of length M, 
The localization energy for each loop is i<i oc ~ — %M whereas the number of loops is (N — m)/M. Thus, the free energy 
of the string of the loops reads 
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where f3% and 0% are constant factors. 

With eqs 01 El El an d El the total free energy (eqQJ reads 
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for strong localization 
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The expression in square brackets of ea 1151 (for the case of the weak localization limit) can be viewed as an inverse effective 
blob size 1/g, i.e. 

1 = 1 _ ai ~- /"R-N 1 ^ 
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The significant difference between the weak and strong localizations which is pointed out in ea 1 151 will become apparent 
in the Sec. 2.2.4. Namely, there we will show that the characteristic translocation time for the weak localization regime 
grows monotone with the block length M, whereas in the strong localization case this dependence passes through the 



maximum. 



2.2 Translocation time 
2.2.1 Preliminaries 

It is well known that the one-dimensional Smoluchovsky problem allows for an exact analytical solution for the mean first 
passage time [24-26]. Let us assume that a stochastic process £(t) which describes the drift - diffusion of a chain in the 
external field U(x) extends between the boundaries x\ or provided that £(0) = ieo- Then the expression for the mean 
passage time t{xq) has the form 

X2 . x' 

//3U(x) f 
dx ' e ~DW) J d y e ~ mv) . 

X KXl 



where D(x) is the corresponding diffusion coefficient. In view of the results of the previous section, we can map the 
problem at hand onto the one dimensional Smoluchovsky problem where the corresponding free energy T{m) plays the 
role of an external potential. The question about the choice of the proper diffusion coefficient D is more open. Considering 
the translocation through a nanopore, Park and Sung [27-29] suggested that D should be identified with the Rouse model 
diffusion coefficient which scales as D — Dq/N. On the other hand, Muthukumar [30] assumed that D is not the diffusion 
coefficient of the whole chain but rather the diffusion coefficient of the monomer which just passes the pore, and hence a 
constant, independent of N (see also an interesting discussion in ref. [31]). We will consider here I? as a free parameter 
and will use a correspondence principle (see below) in order to fix it. We are interested in the average first - passage time 
when the length m changes within the interval 1 < m < N — 1. Using ea 1171 this time can be written in the form 
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We emphasize that this expression is basically valid under conditions when the chain is initially in equilibrium at the 
interface. The completely opposite case is that of a chain that is dragged through the interface while being in the steady 
state conformation typical of a drifting coil. The treatment of the latter case is more complicated and might not be simply 
based on the use of the same free energy unless the dragging motion is pretty slow. 



2.2.2 Kramers' escape time 

One may simplify the double integral of ea 1181 to treat it analytically. To this end only two important terms in ea 1151 can 
be left: the localization energy, —(N — m)/q, and the elastic energy, — (aB /Ty^m^ 1 ^^ 1 : the latter being responsible for 
the detachment. Then the resulting free energy function is shown in Figure 
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Figure 2: Schematic representation of the total free energy (solid line). The dashed line and the dotted line represent the 
localization energy and the elastic energy, respectively. 



As one can see, the resulting barrier closely resembles that of the Kramers' escape problem [24]. In order to roughly 
estimate the escape time and its dependence on B, N and M, we assume that the barrier determines the escape time. 
This means that the slowest part of the process is related to the barrier surmounting, whereas the rest of the translocation 

G 



goes much faster. The free energy functional in this case has the form 

Hm) K _ T E_m T (^y u m «M + i . (19) 

which has a maximum at to* ~ 2i (^gj where the numeric coefficient 71 = 1/(1 + l/v) v . The condition that to* cannot 
be larger than N leads to the restriction (aB/T) > ji/q v N. The corresponding free energy maximum value is given by 

where 72 = 1/0(1 + l/v) 1+v ). 

The integral over to in ea 1181 gets its main contribution at m ~ m*, so that we replace the factor exp{/3J-(m)} 
by exp{#F(m*) - f5\T"(m*)\(m - m*) 2 /2}, where /3T"(m*) ~ -(aB/T) l l u (m*)^/ v ^- x . Now, as to ~ to*, the main 
contribution in the integral over n comes from the area n ~ (we have assumed for simplicity that the low integration 
limit equals zero). Then, exp {(3F(n)} ~ exp{-/3.F(0) - /3.F'(0)ra} where /3T(0) ~ -JV/g and /3.F'(0) ~ 1/<j. After taking 
all these expansions into account the integrals in eq^l can be easily calculated and the result has the following form 
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The following comments are noteworthy in connection with this result. The exponential term has the form of an 
Arrhenius factor exp(TO*/<7) where m* /q is the barrier height (see the solid line in Figure^- The resulting B - dependence 
of hiT appear as lnr ~ T / (aB) which qualitatively agrees with the MC - simulations findings (see below). At B — ► one 
has to* — + N (it can be seen in Figure EJ) and the activation energy is simply equal to the localization free energy of the 
full chain F\ oc ~ N/g (see the corresponding result in ref. [1]) as it should be. 

2.2.3 Unperturbed drift 

As noted above, we will treat the coefficient D in eq[2]as a free parameter which can be fixed by taking advantage of the 
correspondence principle. In our case this principle states that for a sufficiently large external force B the activation barrier 
disappears and one observes a pure drift of the chain. As one can see from ea 1151 this happens when 1/q <C 1. Then the 
interface crossing Td rag is mainly defined by the balance of Stokes friction and external forces, i.e. Co-^ w c.m. = NB, where 
Co is the Stokes friction coefficient and u c .m. is the velocity of the chain's center of mass, given by the formula v c _ m . = -B/Co- 
The effective thickness of the interface in this case is of the order of the blob length, i.e. £ ~ ag v and the time Td rag can 
be estimated as Tdrag — Co Then we can adopt Tdrag as a pre-exponent factor in the Kramers' escape time relation, i.e. 

^mMMm ■ 

In the more general case we should deal with the complete full double integral of equation^] Once again, in this limit 
the second term in ea 1191 prevails and we have 
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It is convenient to substitute the variables in both integrals with: x — (aB/T) 1 ^ m/ 1 ^^ 1 and y = (aB/T) 1 ^ n^ 1 /"^ 1 . 
The resulting expression reads 

a 2 ( T \ 2 /( 1+ ") 7, e~y } J e* 

t= d{^b) j d yjm^j dx ^m^ ■ ( 24 ) 
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By setting the upper limit in the integral over y equal to infinity, we keep in mind that this integral converges as y — > oo 
so that it does not depend on N at large N. This convergence can be easily seen if we take into account the behavior of 
the integral over x for small and large y values, i.e. 
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for y <C 1 



The correspondence to the pure drift regime leads us to the conclusion that the factor in front of the integral in eq| 



should be equal to Td rag (from Tdrag — Co)- This enables us to fix the coefficient D as 
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Finally, we need to point out that the equation for dragging time can be derived differently starting from the first 
passage time formula (II 71) (see Appendix A). 



2.2.4 Numerical analysis 

For a first assessment of our theoretical treatment we present below some numerical calculations based on the Kramers' 
escape time expression (I22fl . First and foremost, we calculate the expression for the blob length g{x,M) which is of 
great importance in ea 1151 for the total free energy. The number of monomers g in the blob can be written as g(x,M) = 
where the crossover selectivity \c has been extensively discussed in ref. [1] and its scaling is 
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Strong localization threshold Xoo 
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Table 1: Weak and strong localization thresholds data 

Xc ~ M ~( 1+ ^/ 2 jV^ 1- ")/ 2 . The substitution of this Xc in 9 leads back to eq^l The characteristic point where the blob 
length becomes of the order of the block length M is the strong localization threshold Xoo ■ Here we will use for Xc and 
Xoo the typical values obtained from Monte Carlo simulations [1] . The corresponding data are given in Table ^ The 
subsequent calculations of the Kramers' escape time in the weak localization regime are done for a selectivity parameter 
X = 1-2. 

We have used these data in the case N = 128 to make a 3D - plot of g which is shown in Figure [3J The contour lines 
in the x ~ -W plane of Figure [21 illustrate the narrow region of M-dependent selectivity strength where a crossover from 
weak to strong localization of the copolymer at the interface occurs. At the critical selectivity % C (M) the blobs reach 




Figure 3: Variation of the blob size g with block length M for various degrees of interface selectivity \. The contour 
lines demark the region of selectivity values in which a gradual transition from weak to strong localization takes place. At 
Xc{M) (the first contour line on the left) one has g ps N (here N = 128) while at \oo (last contour line on the right) the 
blob coincides with the block size M. 



the largest size g » N, which involves nearly all monomers N, whereas with further increase of x the blob size rapidly 
decreases and at the onset of the strong localization, \oo, it reduces to the size of an individual block M. As demonstrated 
below, within the narrow interval of weak localization, Xc{M) < X Xoo(M), confined between the outermost contour 
lines in Figure|3|as well as in the vast region of strong localization, x > Xoo, both our theory and the simulations show a 
qualitatively different behavior of the translocation time t with respect to B and M in the two regions. 

The implementation of g(x, M) in the equation for Kramers' escape time (1221) (where 1/q is given by eq^D) leads to 
a r vs. B relationship which is given in Figure^. One can readily verify that this relationship is of non-Arrhenian type 
which implies that the height of the activation barrier for crossing the interface varies with the variation of the dragging 
force B. Evidently, with growing block size M this non-Arrhenian behavior becomes progressively more pronounced. We 
will see later that this is consistent with MC - results, and that the non-Arrhenian variation of r with B becomes also 
progressively enhanced as the chain length N grows. 

Figure^) displays the variation of r with M for several values of the force B in the regimes of weak and strong 
localization. One can see that the r vs. M relationship differs qualitatively between the two cases. In the weak localization 
regime one finds a very steep (by orders of magnitude) monotonous increase of r with M, i.e. the escape time is strongly 
selective with respect to the block size. In contrast, in the strong selectivity regime, x — 4-0, the change of r with M is 
non-monotonous (it goes through a maximum) and is not so strong. This is a result of the different responses of positive 
and negative loops to the external field B as discussed in Section 2.1.2. 
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Figure 4: (a) t/to vs. dragging force B for several block sizes M. (b) t/tq vs. M in the regime of weak localization, 
X = 1.2, and strong localization, x = 4.0, and several values of the external field B. 



3 Monte Carlo simulations 
3.1 Simulation model 

The off-lattice bead-spring model employed during our study has been previously used for simulations of polymers both in 
the bulk [32,33] and near confining surfaces [34-39]; thus, we describe here the salient features only. Each polymer chain 
contains N effective monomers connected by anharmonic springs described by the finitely extendible nonlinear elastic 
(FENE) potential. 

Ufene = ~y r m [! - R2 ] ■ ( 27 ) 

Here I is the length of an effective bond, which can vary in between £ m i n < £ < £ m axi with £ m in = 0.4, £ max — 1 being the 
unit of length, and has the equilibrium value £q = 0.7, while R = £ rnax —£o = £q — £ m in = 0.3, and the spring constant K is 
taken as K/ksT = 40. The nonbonded interactions between the effective monomers are described by the Morse potential 

U M = £A/{cxp[-2a(r - r mm )] - 2cxp[-a(r - r mm )}} , (28) 

where r is the distance between the beads, and the parameters are chosen as r min = 0.8, €m = 1, and ol — 24. Owing 

to the large value of the latter constant, Um{t) decays to zero very rapidly for r > r m i n , and is completely negligible for 

distances larger than unity. This choice of parameters is useful from a computational point of view, since it allows the use 

of a very efficient link-cell algorithm [40]. From a physical point of view, these potentials eqs [23 1281 make sense when one 

interprets the effective bonds as Kuhn segments, comprising a number of chemical monomers along the chain, and thus 

the length unit £ max — 1 corresponds physically rather to 1 nm than to the length of a covalent C — C bond (which would 
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only be about 1.5 A). Since in the present study we are concerned with the localization of a copolymer in good solvent 
conditions, in cq[5S|we retain the repulsive branch of the Morse potential only by setting Um{t) = for r > r m i n and 
shifting U M (r) up by e M - 

The interface potential is taken simply as a step function with amplitude \i 

U mt (n,z)=l (29) 
[<7(n)*/2, z<0 

where the interface plane is fixed at z — 0, and, as explained before, o~(n) = ±1 denotes a "spin" variable which distinguishes 
between P- and H- monomers. The energy gain of each chain segment is thus —\, provided it stays in its preferred solvent. 

During each Monte Carlo update, a monomer is chosen at random and one attempts to displace it randomly by 
displacements Ax, Ay, Az chosen uniformly from the intervals —0.5 < Ax, Ay, Az < 0.5. The transition probability 
for such an attempted move is simply calculated from the total change AE of the potential energies defined in eqs |27| - 
1291 as W = exp(— AE/ksT). According to the standard Metropolis algorithm, the attempted move is accepted only if W 
exceeds a random number uniformly distributed between zero and unity. Since our potentials are constructed such that 
chains cannot intersect themselves in the course of random displacement of beads, one does not need to check separately 
for entanglement restrictions. 

One should point out that the external field B, which imposes a bias in the random hops performed by the chain 
monomers, should be taken in the simulation sufficiently weak [41], BAZ <C feT. Indeed, the field B whose only 
component is directed along the negative Z-axis enters as an additional term in the Boltzmann factor, SE = BAZ, which 
makes the energy a steadily decreasing function of Z. The Metropolis algorithm thus drives the polymer chain towards an 
unreached (and unreachable) minimum which produces a uniform drift downwards. The average jump distance in a dilute 
system (AZ) w 1/4, no matter how strong the applied field is chosen [41]. Thus, for large B the probability of a jump 
along the field quickly saturates at unity so that any further increase of B would not really lead to higher drift velocity of 
the polymer. 

In the course of the simulation the starting configuration of the copolymer is relaxed for a period of time and then 
placed at the ceiling of our simulation box before the field B acting in vertical direction is switched on. The onset of the 
translocation process is defined as the time when the first monomer of the drifting chain touches the interface and the 
moment when the last monomer leaves the upper half of the box then denotes the end of the process. Since the time r 
during which the chain remains in contact with the interface fluctuates considerably, for all chains of length 16 < N < 256 
we take averages over 1000 simulation runs. 

3.2 Discussion of Monte Carlo results 

In section [5] we assumed a mechanism of copolymer detachment from the interface which have lead us to predict a non- 
Arrhenian variation of the time r with the imposed field B. We demonstrated that a selective liquid- liquid interface can 
be very sensitive with respect to the block size M of a multiblock copolymer chain. On the other hand, as shown in 
Appendix^ the N dependence of r is rather weak and nearly vanishes in the limit of long chains. In order to test the 
predictions of the analytic theory and deepen our understanding of the overall picture of translocation through a liquid 
interface, we report in the present section the results of sonic extensive Monte Carlo simulation by means of the model 
described in Section ^ 
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Figure 5: Master plot of the distribution of capture times P(t) vs. t/(t) for a copolymer chain with N = 128, M = 4 
at x = 1-5 for three different fields -B, as indicated. The time intervals have been scaled by the value of (r) appropriate 
for the respective value of the field B. A full line denotes the best fit of a r-distribution. 

In FigureElwe show the distribution of capture times P{t), measured for a copolymer in the weak localization regime, 
for three different field strengths. Despite the considerable scatter in data which remains after 1000 measurements, 
evidently these distributions are shown to collapse on a single "master" curve when scaled by the mean capture time 
(t(B)). The solid line in Figure|S]corresponds to a r-distribution. This distribution arises naturally in processes for which 
the waiting times between Poisson distributed events are relevant. 

One of the main results of our computer experiments is displayed in Figure El where the variation of (r) with the bias 
B is plotted for chains of different length N and fixed block size M = 4 in semi-log coordinates. A general feature, evident 
from FigureEKj is the well expressed non-Arrhenian (r) vs. Bl/ksT relationship which confirms our expectations from 
Section [5]- see Figure 0] >From Figure it becomes evident that this non-Arrhenian type relationship crosses over to 
Arrhenian one for shorter chain (N = 16 , 32) where insufficient length prevents the formation of (large) blobs of size g as 
expected at criticality Xc — 1.5. For sufficiently long copolymers, N > 64, in contrast, the capture time r depends only 
weakly on chain length N and the t vs. B relationship remains qualitatively unchanged. This suggests that the activation 
barrier, which has to be overcome when a chain is detached from the interface by the field, grows with decreasing drift 
velocity. One could imagine that in a strong field B the copolymer moves so fast that the translocation time is much 
shorter than the time it takes for the polymer to localize at the phase boundary and attain its lowest free energy. For 
small values of Bl/ksT, in contrast, the copolymer has sufficient time to attain full equilibrium at the selective interface, 
and the corresponding potential well in which the chain resides gets deeper. While the non-Arrhenian r vs. Bl/ksT 
relationship is also observed in the case of strongly selective phase boundary, Figure 03, qualitatively the variation of r 
with chain length N looks quite different. While at both ends of the interval of variation of B the capture time practically 
does not depend on N, in the intermediate region, 0.4 < BL/ksT < 0.8, the longer chains spend less time at the interface 
than the shorter ones! At present we lack a reliable explanation for this effect. One might assume, however, that in this 
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Figure 6: (a) Variation of r vs. dragging force B for chain lengths N = 16 , 32 , 64 , 128 , and 256 in the weak localization 
regime \ = 1-5- (b) The same as in (a) in the regime of strong localization, \ = 4-0, 



interval of intermediate B intensities the process of copolymer detachment is governed by the probability for the formation 
of a critical protuberance. Such a probability should then be proportional to the total chain length N. Considerations of 
this kind indicate, however, the necessity of a more detailed information about the underlying mechanism of copolymer 
detachment and its temporal evolution. Therefore, in Figure we present a series of snapshots which capture the time 
history of detachment by depicting the instantaneous positions of all repeat units of the chain at four successive times after 
the onset of one particular translocation. One can clearly see in Figure0that after an initial latent period a protuberance 
is formed at the free (dangling) end of the chain between monomers 120 < n < 128 which from the start is immersed in its 
preferred solvent. After 20000 MCS this end evolves into a progressively growing protuberance which eventually tears off 
the copolymer from the interface and drags it into the lower half of the simulation box. Our observation show, however, 
that such protuberances form occasionally everywhere along the backbone of the chain and then die out with time or 
grow spontaneously like an avalanche. Thus the analogy with the mechanism of detachment, envisaged in our theoretical 
treatment, Section |2 should be apparent. We note here that in order to check the role which the dangling end of the 
chain plays in assisting the detachment process we also made test runs with ring copolymers (see below) under otherwise 
equivalent conditions and found no qualitatively different behavior when dangling ends are eliminated. 

An important test for the theory, as developed in Section [21 is represented by the analysis of the block size dependence 
of the translocation time r(M). In Figure |SJi we demonstrate the extremely strong sensitivity of t on M for the case of 
weak selectivity, x = 1-5, which leads to a steady increase in r by more than a decade with doubling the block size. Figure 
IHt also suggests that this sensitivity gets fully pronounced for sufficiently long chains, N > 64, which we attribute to the 
formation of sufficiently large blobs unrestricted by the finite size of the copolymer. The variation of the translocation 
time t with block size M for a fixed field strength B = 0.3 ,0.5 ,0.7 is shown in|HjD where the steep (nearly exponential) 
increase is evident. Qualitatively the observed behavior is well in line with the one predicted by our theory and shown in 
Figure OJl. 
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Figure 7: Snapshots showing the time evolution of the copolymer detachment from a liquid-liquid interface for a chain 
with TV = 128 and M — 8 (different blocks are represented by solid and dashed lines) at \ = 1-5- Displayed are the Z - 
coordinates of all 128 monomers with respect to the position of the phase boundary at Z = 0. The times shown are at 
325 (a), 20325 (b), 100325 (c), and 200325 MCS (d) after the start of the translocation process. 

For the case of strong selectivity, \ = 4.0, again we find a qualitatively different r vs. B behavior, see Figure |31 which 
agrees quite well, however, with the theoretical predictions - Figure 0Jd. For a given intensity of the dragging force one 
observes a non-monotonic change of r with increasing block size M with a local maximum around M k, 4. Thus one could 
conclude that our computer experiments provide a strong support for the theoretical model, developed in Section [21 on the 
ground of simple scaling considerations. 

Eventually, before the end of this section we would like to point out several additional results which are believed to 
underline the consistency of the present combination of analytic and simulational results. An important question concerns 
the impact of the starting conformation of the copolymer on the observed non-Arrhenian t vs. Bl/ksT relationship. To 
this end in Figure HHI we display a comparison of the measured t(B) dependence for chains which are localized and fully 
equilibrated at the liquid-liquid interface before the drag force is switched on, and chains with identical N and M which 
traverse the phase boundary "on the fly". As one may readily verify from Figure ITU1 for strong fields, B > 0.03, the chains 
which start from equilibrated position at the interface need somewhat shorter time r to cross it than those arriving from 
the upper half of the simulation box. This appears reasonable since the former chains "save" the time interval between 
first touch of the interface whereby their center of mass still lags behind at Z cm f« R g above the separation line, and the 
moment when their center of mass arrives at the interface, Z cm « 0. For sufficiently small forces, however, one can see 
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Figure 8: (a) Variation of the mean translocation time r with field strength B for copolymers of length N — 128 and 
different block size M in the weak localization regime \ — 1-5. (b) Variation of r with block size M for B — 0.03, B = 0.05 
and B = 0.07. An exponential fit is indicated by a dashed line. 

in Figure ITU1 that both curves merge, that is, the time of residence at the interface is so long that the small temporal 
advantage for a chain starting from the interface is completely lost. Both in the case of "equilibrated" and "on the fly" 
chains, however, qualitatively the r vs. Bl/ksT relationship and its non-Arrhenian type remains unaltered which justifies 
the use of equilibrium free energy expressions for the copolymer in our theoretical treatment. 

Another question which pertains to the specific mechanism of copolymer detachment from the interface concerns the 
role of free (dangling) ends which are expected to facilitate the process of detachment. A comparison between linear 
chains of a given length, say N — 128, and ring copolymers of the same length might thus expose the role of free chain 
ends. From Figure ITUI one can see that for strong forces, B > 0.025, the translocation time is about 50% shorter in the 
case of ring whereas for B < 0.025 it becomes gradually considerably larger than for linear chains with dangling ends. 
Indeed, for strong fields when the objects do not have sufficient time to localize at the interface the ring polymer crosses 
faster the border line between the solvents because its size is smaller than that of the linear chain. In this regime the 
polymers move as almost compact objects and the presence of free ends does not affect the the mechanism of interface 
crossing. In contrast, for weak forces the polymers have sufficient time to localize and equilibrate on the interface while 
being dragged through it. The effort to tear a dangling end off the interface would then be only half of what is needed 
to create a protrusion in a ring copolymer since the protrusion is hinged by two ends to the interface. Thus it should 
not be surprising that the detachment of ring copolymers take considerably longer than the that of the linear ones, as 
demonstrated by Figure E3 

A special case which has not been in the focus of the present investigations albeit it may serve as a check for our 
dynamic Monte Carlo simulation is shown in the inset of Figure ITUI There we present the variation of translocation time 
t with chain length N for a telechelic copolymer with fixed block size M = 4 whereby a constant drag force B = 0.9 is 
applied only to the head monomer of the chain. Evidently, one recovers a quadratic increase of the translocation time with 
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Figure 9: (a) Variation of the mean translocation time r with field strength B for copolymers of length N — 128 and 
different block size M in the strong localization regime x = 4.0. (b) Variation of t with block size M for B = 0.3 , 0.5 , 0.7. 
As expected, cf. Figure 0Jd, the r vs. M relationship is for strong selectivity not a steady one. 



length N, r oc N 2 , which is qualitatively different from our data when the external field acts on all monomers of the chain. 
One can easily understand this observation, however, noting that for a telechelic chain the friction force grows linearly 
with the number of beads N, meaning that the drift velocity slows down as N^ 1 . In addition, the characteristic size of 
a stretched chain should be proportional to N too. The combination of both factors should yield a quadratic growth of 
r with N, and this is indeed nicely confirmed by the simulation. This last result gives us an additional confidence in the 
consistency of our simulation model and the observed results. 

4 Conclusions 

In the present work we have studied the transition of field-driven regular block copolymers through a selective liquid-liquid 

interface. We have developed a simple theory based on scaling considerations which appears capable of catching the main 

features of the phenomenon, namely, the pronounced non-Arrhenian dependence of the mean translocation time r on field 

intensity B, and its extreme sensitivity with respect to the copolymer architecture, that is, the block size M. It is based 

on the notion that the mechanism of copolymer detachment from the interface can be modeled as an activation process 

whereby the copolymer is pulled out of a potential well, determined by the free energy of a localized equilibrium chain, 

by a progressively growing protuberance. The height of the activation barrier is shown to depend dramatically on the 

field strength B - for large values of Bl/ksT the barrier height disappears and the chain is dragged without resistance 

through the selective phase boundary whereas for weak fields the free energy of compressed and stretched blobs determines 

r. The mean translocation time itself is tackled as a mean first passage time using the Kramers approximation and thus 

predictions for qualitatively different behavior of r in the regimes of weak and strong selectivity are derived. 

In order to check the analytic theory and provide a more comprehensive understanding of the process under investi- 
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Figure 10: Variation of r with field strength B for a chain with length N = 128 and block size M — 4. Circles denote 
a ring copolymer while squares mark linear chains which start from equilibrated initial conformations localized at the 
interface (full symbols), or enter the interface "on the fly" (empty symbols). The inset shows the variation of r with N for 
a telechelic chain where the field B = 0.9 is applied only at the head monomer of the chain. The straight line denotes a 
power law, r oc N@ , with an exponent (3 = 2.02 ± 0.04. 



gation, we have reported on extensive simulations carried out by means of a dynamic Monte Carlo model which strongly 
support our findings. The simulations confirm the characteristic non-Arrhenian type of r vs. Bl/ksT relationship while 
revealing distinct differences for the cases of weak and strong selectivity of the liquid-liquid interface. One remarkable 
feature which is reproduced by the computer experiments is the strong sensitivity of the translocation time with respect to 
the block size M whereas the total chain length N turns to have little effect on r. This sensitivity suggests the interesting 
possibility to use selective liquid-liquid interfaces as a new type of chromatographic tool whereby one can "sieve", that is, 
separate and analyze, complex mixtures of copolymers with respect to block size M. Recently a novel method to charac- 
terize individual blocks by means of Liquid Chromatography at the Critical Condition (LCCC) [42] has been proposed. 
In the present work we discuss another possibility of a new type of chromatography, based on the afore mentioned results. 

One should note that in the rich behavior reported above we have not included results pertaining to random copolymers 
at penetrable selective interfaces where the range of sequence correlations plays a role similar to that of the block size 
M. Moreover, here we have totally ignored the possible effects of hydrodynamic interactions on the translocation kinetics. 
One might expect that the existence of a sharp interface, separating the two solvents, as well as the probable significant 
difference in the solvents viscosities would generally suppress hydrodynamic effects across the phase boundary. Of course, 
the problem becomes more complex if the presence of a copolymer at the interface changes the interface properties. These 
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are important complications which deserve a thought investigation. Nevertheless, we could expect that the main finding 
about the strong sensitivity of r with respect to M is unchanged. Eventually, this effect is evident from the localization 
free energy of the string of blobs (see eas !12l - [T^|l . where the main term grows as -Faring ~ —TNAI 4 , i.e. strongly depends 
on M. Clearly, the verification of these predictions requires considerable computational efforts and therefore is on the 
agenda in our future work. 

One should also keep in mind that all results have been derived and checked against Monte Carlo simulations within 
the simplest model of an interface of zero thickness. It is conceivable that by allowing for the presence of an intrinsic 
width of the interface as well as for a more realistic description involving capillary waves, etc., many additional facets of 
the general picture should become clearer. While these and other aspects imply the need of further investigations, i. e., 
by means of single-chain laboratory experiments, we still believe that the present study constitutes a first step into a 
fascinating world which might offer broad perspectives for applications and development. 
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A Calculation of the pure drag time 



In this Appendix we calculate the characteristic drag time for a polymer of length N pulled by an external force B along the 
z - axis. The calculation is based on the first passage time formula l|17|) which enables also to estimate 1/N - corrections. 
In this case the energy field in ea 1171 has the form 
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(30) 



The polymer chain center of mass starts from z = at t — 0; we want to estimate the value of the time at which the 
center of mass of the chain passes z = L. According to eqElthe average time of the first passage reads 
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where we have used for the diffusion coefficient its Rouse expression D — T/C,qN. The straightforward calculation of the 
integrals in ea 1311 yields 
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One can see that for a large chain , i.e. at N ^> 1 and for a characteristic length of the order of the blob length, i.e. at 



L ~ ag v we go back to the simple expression raraj 
factors proportional to powers of [1 — Tj (BLN)]. 



Co- The finite chain size correction can be expressed in terms of 
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